Zero-bias anomaly and Kondo-assisted quasi-ballistic 2D transport 
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Nonequilibrium transport measurements in mesoscopic quasi-ballistic 2D electron systems show 
an enhancement in the differential conductance around the Fermi energy. At very low temperatures, 
such a zero-bias anomaly splits, leading to a suppression of linear transport at low energies. We also 
observed a scaling of the nonequilibrium characteristics at low energies which resembles electron 
scattering by two-state systems, addressed in the framework of two-channel Kondo model. Detailed 
sample-to-sample reproducibility indicates an intrinsic phenomenon in unconfined 2D systems in 
the low electron-density regime. 
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Zero-bias anomaly (ZBA) in the non-equilibrium char- 
acteristics of ballistic systems provides additional insight 
into the coupling mechanism of a conduction electron 
with its surroundings. For example, a maximum in the 
differential conductance (dl/dVsu) at the Fermi energy 
(E-p) in case of tunnelling via magnetic impurities 
or through artificially confined quantum dots 0, 0] , has 
been explained by a Kondo-like antiferromagnetic cou- 
pling of the electron to localized impurity spin. Similar 
enhancement observed in clean quantum point contacts 
in semiconductor heterostructures has led to much con- 
troversy regarding the spin-structure and effects of lateral 
confinement in such systems Q. On the other hand, in 
quasi-ballistic metallic nanobridges the ZBA shows in a 
local cusp-like minimum in dl/dVsu at Ep, which has 
been interpreted in terms of nonmagnetic two-channel 
Kondo (2CK) framework arising from the interaction of 
the electrons with local two-state atomic defects . The 
interest in this lies in the prediction that in particle-hole 
symmetric case, such a model flows to a T = fixed 
point, giving rise to quantum critical behavior and non- 
Fermi liquid effects ja, . 

Investigation of nonequilibrium ballistic transport in 
spatially extended 2D systems is primarily impeded by 
greater scattering from background disorder. The small 
level spacing requires relatively larger sample dimensions, 
resulting in significant momentum relaxation during the 
course of transport. Even though recent nonequilibrium 
studies in mesoscopic 2D electron systems (2DES) have 
shown an unexpected ZBA at low background disor- 
der 8] , the mechanism of transport in this regime remains 
unclear. Here we report the experimental observation of a 
new and unexpected feature in the nonequilibrium quasi- 
ballistic 2D transport. At zero magnetic field in uncon- 
fined samples of mesoscopic dimensions, we find a strong 
zero-bias enhancement in dl/dVsu, which splits at very 
low T (< 150 mK) by a gate voltage (V^)-dependent mag- 
nitude. This modifies the linear transport properties by 
introducing intermittent low-energy nonmonotonicity in 
the linear conductance G (= dl/dVsn at Vsd = 0, where 
Vsd is the source-drain bias), as a function of both T and 



in-plane magnetic field {Bn). At low Vsd, an intriguing 
scaling behavior of dl/dVsT) in T is also observed, which 
indicates a 2CK-type scattering of electrons in these sys- 
tems at low T. 

The devices were fabricated from Si i5-doped 
GaAs/AlGaAs heterostructures. In order to minimize 
the disorder arising from Coulomb potential of ionized 
dopants, we used a thick layer (« 80 nm) of undoped Al- 
GaAs spacer, and adopted a slow cooling procedure || 
to maximize the correlations in the donor layer. The as- 
grown mobility of the wafers range over /i~l-3x 10 6 
cm 2 /V s, depending on the donor (ns) and electron (n s ) 
density. Over the experimental temperature range, the 
momentum relaxation rate was typically r _1 < (0.1 — l)x 
ksT/fL indicating a quasi-ballistic nature of the trans- 
port |9j. The elastic scattering length A ~ vpr ~ 10 /im, 
where i>f is the Fermi velocity, provided an upper cut- 
off to the device dimensions, restricting the number of 
elastic scattering events to very few or none. Measure- 
ment of dl/dVsT) was carried out using standard mixed 
low-frequency ac/dc technique with ac excitation <gC k^T 
at each T inside a top-loading dilution refrigerator with 
base electron temperature 32 mK. 

Fig. 1(a), (b) and (c) show the B\\ = nonequilibrium 
characteristics of three samples with different dimensions 
and doping profile. We varied both the total geometrical 
area, as well as the aspect ratio in samples T45 (2 x 8/im 2 , 
n s 2.5 x 10 12 cm- 2 ), T46 (5 x 8^m 2 , n& » 1.9 x 10 12 
cm~ 2 ) and T48 (5 x 5^m 2 , n s » 0.9 x 10 12 cm~ 2 ). The 
dimensions are chosen in order to ensure small single- 
particle level spacing with Ae x ,Ae y ~ h 2 /8m*L 2 <C ksT 
for each T. While no consistent feature in G as a function 
of V g could be identified (not shown), the ZBA around 
Vsd = is evident for all samples . We focus on some 
of the common features of the ZBA shown in Fig. 1: (1) 
At | Vsd | < 0.15 — 0.2 meV, dl/dVsD shows an enhance- 
ment at all V s . This energy scale varies only weakly with 
V^, and also -C Ep at all n s . (2) Over certain ranges of V g , 
the ZBA splits by 2A, leading to a double-peak structure 
around E-p. (3) A is oscillatory in V s , and when normal- 
ized for series conductance, in particular for G > 2e 2 //i, 
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(a) Sample: T45 



Sample: T46 



(c) Sample: T48 
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FIG. 1: (Color online) The zero-bias anomaly at T ~ 32 mK 
and Bu =0: (a)-(c) Nonequilibrium characteristics in three 
different samples. Each differential conductance trace is ob- 
tained for a fixed V g . (See text for sample details) (d),(e) 
B| | -dependence of the nonequilibrium characteristic at the 
single peak and double peak regime (T = 32 mK). Traces 
are vertically shifted for clarity. The vertical bars denote a 
conductance change of 0.2 x e 2 /h. 




1.0 



g 0.8 

T3 



-100 -50 50 -100 


-50 50 


eV s A T 


eV S D /k B T 


(e) ° 00 °^ 




9 peak 1 




O peak 2 





0.01 



0.1 



T(K) 



the splitting is strongest when G lies close to an even 
integral multiple of e 2 /h, becoming unresolvable around 
odd multiples. (4) While no clear dependence of the en- 
ergy scales on size/shape of the samples was observed, 
the amplitude of the ZBA was found to depend on n s . 
Note that the ZBA is strongest in T48 (n s = 6.5 x 10 10 
cur 2 at Vg = 0), and weakest in T45 (ra B = 10.1 x 10 10 
cm~ 2 at V" g =0). Such weakening is also observed in a 
given sample as Vg (or n s ) is increased. 

The single- and double-peak regions show distinct be- 
havior in the presence of external in-plane magnetic field 
(_B||). As shown in Fig. 1(d) and (e), while the single- 
peak splits by the Zeeman energy Az = 2g*^B-E>|| from 
_B|| =0, Zeeman splitting of the double peak appears only 
at relatively large (sample-dependent) i?||-scale (~ 0.5 — 2 
Tesla) 8J. Such a splitting of ZBA is taken as a distinc- 
tive feature of Kondo-type dynamics Q, Q- The effec- 
tive g-factor g* was found to be both sample- and rig- 
dependent g* /gb ~ 1 — 3, where \gb\ = 0.44, consistent 
with exchange-induced enhancement at low n s 0, flll | . 

Kondo-type dynamics also lead to a suppression of the 
ZBA with increasing T. Indeed, such a suppression is 
observed in all our samples as shown in Fig. 2a for nine 
different V^'s in T46. In case of single peaks (sets 2, 3, 
7 and 8), the suppression is monotonic as T is increased 
from « 30 mK to » 800 mK. For A ^ 0, the T depen- 
dence of dl/dVsT) is nonmonotonic in T close to Ep (see 
sets 1, 4, 5, 6, and weakly in 9). Fig. 2b shows the T- 



FIG. 2: (Color online) (a) Suppression of the zero-bias 
anomaly with increasing temperature (T). Each of the nine 
sets was recorded at a fixed gate voltage (Vg) and Bu — 0. 
(b) Dependence of the linear conductance G on T at different 
Vg. The length of the vertical bar denote 0.1 x e 2 /h. (c),(d) 
Collapse of AG onto a single trace at Iow-Vsd as a function of 
eVsD /ksT, at two regimes of zero-bias conductance (G). The 
insets show the unsealed original data, (e) T-dependence of 
individual peaks in a double-peak pattern. 

dependence of G from the traces in Fig. 2a. Qualitatively, 
the high-T regime (typically T > 300 mK) is similar for 
all cases, where G shows a "metal" -like decrease with in- 
creasing T. For T < 100— 150 mK, the increasing behav- 
ior of G with T appears intermittently, at values of V g 
where dl/dVsu shows resolvable double-peak structure 
(Fig. 2b and 3b-e). This reentrant nature is crucial, be- 
cause even if the high-T suppression of G is attributed to 
various combinations of phonon contribution, the inter- 
action correction, or the Altshuler-Aronov-type correc- 
tion from electron scattering by Friedel oscillations [l^ . 
the repeated change in the sign of dG/dT is clearly in- 
consistent with standard weak-localization/interaction- 
based mechanisms having a single transition from local- 
ized to "metallic" state transport. However, in the Zee- 
man split regime at high Bu, G increases with T at all 
Vg, indicating standard localized state transport. 

Most striking aspect of the double-peak structure is 
the collapse of AG(V SD , V S ,T) = dI/dV SD (V S B,V gl T) - 
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G(Vg,T), onto a single Vg-dependent trace when Vsd is 
scaled by T. In Fig. 2(c) and (d) this is illustrated for 
ZBA's at two different G (and hence V g ), indicating a 
common underlying mechanism. The insets contain the 
actual data prior to substraction of G(V g ,T), showing 
the thermal broadening of dl/dVsn, and the suppression 
of the ZBA splitting with increasing T. In Fig. 2e we 
show the T-dependence of individual peaks in dl/dVsu 
in the double-peak region at the given V g . For compari- 
son, the individual T-dependences are normalized by the 
dl/dVsD at lowest T. Note the approximately logarith- 
mic T-dependence of both peaks at T > 150 mK, which 
is also expected in the Kondo-framework. 

The low- Vsd scaling of dl/dVsu indicates nonequilib- 
rium conductance of the form, 



dl 
dVs^ 



(Vsb,V s ,T) = G(V S ,T) + AT<*T 



V*, 



eVsp 
k B T 



(1) 



where A is a phenomenological constant, and a — gives 
the best collapse onto the single sample-dependent func- 
tion F (Vg, eVsT> / k^T) . The dependence of T on V g will 
be discussed later. The scaling form of Eq. has been 
observed in the context of the scattering of electrons 
from bistable systems 0, and addressed in a two- 
channel Kondo (2CK) framework 0,0- In quasi-ballistic 
nanostructures, 2CK dynamics may arise from two pos- 
sibilities: (1) Electron scattering off systems with quasi- 
degenerate orbital states acting as pseudospins, while the 
real spins act as the channel index. The nature of this 
orbital degeneracy can however vary from equivalent lat- 
tice defects in metallic nanobridges , to singlet-triplet 
degeneracy in multilevel quantum dots [Til Hq. (2) An 
underscreened high-spin (S) gro und state coupled to M 
(< 2S) conduction channels |lq |. When parametrized in 
terms of a two-impurity Kondo problem, two spins Si 
and 52 (S — S\ + S2) interact with conduction electrons 
with antiferromagnetic exchange parameters J\ and J2 
{J2 7^ Ji), and a direct exchange I (J\, J 2 S> /) [Tjj. A 
two-stage screening process can then decouple the spins 
from the conduction band by forming S1-S2 singlet if I 
exceeds some critical magnitude. Since both stages can 
be suppressed by lifting the spin degeneracy with Zeeman 
energy, a distinctive feature of this mechanism is the non- 
monotonicity of G in both T and B\ | , with a quantitative 
correspondence between the respective energy scales |l6j . 

In Fig. 3, we compare the T-dependence of G to its In- 
dependence in T45 (on a separate cooldown). Four rep- 
resentative V g , with corresponding noncquilibrium traces 
at T ss 32 mK and Btt = T, are identified as V g i to V g 4 
in Fig. 3a. Figs. 3b-e show the (B\\ = 0) T-dependence of 
G at these V g s, while Figs. 3f-i show the B\\ -dependence 
at the base T w 32 mK. We note that apart from the 
qualitative nonmonotonic behavior of G as a function of 
both T and B\ \ at V g \ and V g s , there is also a quantitative 
agreement in the energy-scales over which the double- 
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FIG. 3: (Color online) Correspondence between the source- 
drain bias, thermal and Zeeman energy scales: (a) Nonequi- 
librium traces in T45 identifying two single and two double- 
peak regions at gate voltages V^,i to V S 4. (b)-(e) Temperature 
(T)-dependence of the linear conductance G. (f)-(i) In-plane 
magnetic field (-B||) dependence of G at corresponding V g s. 



peak structures are suppressed. For example, as the half 
width A/2 reduces from « 0.03 meV at V gl to « 0.017 
mV at V g 2 , we find a corresponding decrease in the char- 
acteristic thermal and Zeeman energy scales (denoted by 
the vertical arrows) from ~ 0.023 meV and ~ 0.028 meV 
respectively at V s \, to ~ 0.011 meV and ~ 0.02 meV 
respectively at V g 2- Since B|| is applied in the plane of 
the 2DES (to an accuracy better than 0.2°), conventional 
weak-localization effects are unlikely to contribute to the 
observed behavior [19j . Coupling of B\ | to orbital degree 
of freedom through spin-orbit interaction [2(| or finite 
thickness effect |21| are also excluded since neither antilo- 
calization to weak-localization crossover, nor a stronger 
suppression of G with increasing B^ at lower n s were 
observed. 

The Kondo-type enhancement in mesoscopic 2D non- 
linear transport in the ultra-clean limit indicates an unex- 
pected scattering mechanism of the lead electrons. More- 
over, intermittent splitting of the ZBA implies a two- 
state nature of the scatterer that becomes resolvable 
only at very low T (< 150 mK). In view of the quasi- 
ballistic nature of transport, the non-equilibrium con- 
ductance can represent the tunnelling characteristics be- 
tween the source and drain across the potential barrier 
formed by the gate. Presence of localized magnetic states 
within the barrier region of traditional tunnel junctions 
has been shown to result in Kondo-type enhancement in 
the tunnelling conductance 0. In nanostructures fabri- 
cated from MBE-grown high-quality GaAs/AlGaAs het- 
erostructure, localized acceptor sites close to the sys- 
tem can significantly modify the transport through ca- 
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FIG. 4: (Color online) Scaling of the low-energy nonequi- 
librium characteristics with a gate-voltage (V^)-dependent 
energy-scale F(V g ) at T » 32 mK for (a) T46 and (b) T48. 
Insets show the unsealed, but zero-shifted, traces. Both sets 
were obtained for G around 4e 2 //i in Fig. 1. (c),(d) Variation 
of F(V^) as a function linear conductance G in the T46 and 
T48 respectively. Solid lines are guide to the eye. 

pacitive coupling |22j. This could lead to a Kondo- 
assisted tunnelling, similar to that observed in metallic 
nanobridges H, and also in semiconductor-metal junc- 
tions 24]. Arguments against such a scenario are, (1) 
reproducibility, and insensitivity of the relevant energy 
scales to sample-specific details, and (2) enhancement in 
the amplitude of ZBA at lower disorder. 

Alternatively, dependence of the effect on n s , i.e. rel- 
atively strong amplitude of ZBA in samples with low-n s , 
and also weakening of the effect with increasing n s within 
a given sample, indicates a possible role of Coulomb in- 
teraction. Typically in our systems, the low magnitude 
of n a (~ 0.8 — 3 x 10 10 cm -2 ) results in a large inter- 
action parameter, r s — 1/ a^,^/i:n s ~ 3.5 — 6, where <xb 
is the effective Bohr radius. This corresponds to an ex- 
change energy that can lead to strong many-body spin 
fluctuations with magnetic moment 3> 1/2 |25j. Under- 
screening of such spin fluctuations when Kondo-couplcd 
to lead electrons can result in split Kondo resonance, 
which would be strongly V^-dependent. Exchange split- 
ting of Kondo-resonance in presence of itinerant electron 
ferromagnetism has been observed experimentally |2fij . 
Here, the suppression of the double peak pattern in Bu 
provides further evidence of the dynamics to be related 
to real spin, rather than pseudospin resulting from quasi- 
degenerate orbital states. Indeed, similar ZBA observed 
in quasi- ID quantum point contacts have been inter- 
preted in terms of a dynamic spin-polarization of the 
electrons within the ID channel |4|. 

Finally, we discuss the intermittent collapse of A ob- 
served in Fig. 1. If the scaling relation of Eq. 2] originates 
from proximity of a 2CK fixed point H 03, the low-T 
VJr-dependence of A indicates that such scaling should 



also be possible in terms of a V^-dependent energy-scale 
r(V^). This implies an asymptotic form of the scaling 
function T -> T[V su /T(y e )] as T -► 0. Theoretically, 
r is analogous to the energy asymmetry in case of or- 
bital degeneracy, or the direct exchange parameter / in 
the two-impurity Kondo parametrization [ljj [lfj, UJl UM > 
which depends implicitly on V g through the exchange in- 
teraction magnitude J j27}- Indeed, in Fig. 4a and 4b 
we have shown the low-energy scaling of AG(Vsd, ^g) = 
dI/dV S D(VsD,V e ) - G(Vg) in T46 and T48 at the base 
T 32 mK. In the absence of an absolute scale, we 
have expressed T(V g ) with respect to its minimum value 
(r(V^) mm ) observed around G ~ 3.5 x e 2 /h in both sam- 
ples (Fig. 4c and 4d). The scalability of AG close to 
the collapse of A strongly indicates the possibility of a 
quantum critical dynamics. However, a satisfactory ex- 
planation of its reentrant nature and the fundamental 
mechanism of the 2CK-type scattering forms the basis of 
ongoing investigations. 
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